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Abstract 

We calculate the decay rates of B mesons into P-wave charmonium states using 
new factorization formulas that are valid to leading order in the relative velocity of 
the charmed quark and antiquark and to all orders in the running coupling constant 
of QCD. We express the production rates for all four P states in terms of two non- 
perturbative parameters, the derivative of the wavefunction at the origin and another 
parameter related to the probability for a charmed-quark-antiquark pair in a color-octet 
S-wave state to radiate a soft gluon and form a P-wave bound state. Using existing 
data on B meson decays into Xd to estimate the color-octet parameter, we find that 
the color-octet mechanism may account for a significant fraction of the Xd production 
rate and that B mesons should decay into Xc2 at a similar rate. 



The production rate of quarkonium states in various high energy physics processes 
can provide valuable insight not only into the interactions between a heavy quark and anti- 
quark, but also into the elementary processes that produce the QQ pair. The spin-1 S-wave 
resonances, like the J / ip of charmonium, are of special experimental significance because they 
have extremely clean signatures through their leptonic decay modes. Because a significant 
fraction of the ip's come from the decays of the P-wave Xcj states, an understanding of the 
production of P-wave resonances is necessary in order to understand inclusive ip production. 
The P states are also important in their own right because they probe a qualitatively dif- 
ferent aspect of the QQ production process. While the S states probe only the production 
at short distances of a QQ pair in a color-singlet state, the P states, as we shall show, also 
probe the production of a QQ pair in a color-octet state. 

One of the simplest production processes for charmonium states is the decay of a B 
meson or baryon. For B~, B°, B s , and A& (but not for B c ), the c and the c that form the 
charmonium bound state must both be produced by the decay or annihilation of the b quark. 
Since this is a short distance process which occurs on a length scale of order 1/M&, where Mj, 
is the mass of the b quark, it should be possible to apply perturbative QCD to calculate the 
inclusive decay rate into a particular charmonium state. If we neglect contributions that are 
suppressed by powers of A/M&, where A is a typical momentum scale for light quarks, the 
decay rate of a B hadron is given by the decay rate of the b quark, with the light antiquark 
in the B meson and the light quarks in the B baryon treated as noninteracting spectators. 
In the case of the hadronic and the semileptonic decay rates, the leading corrections are 
suppressed by two powers of A/M& (Ref. [I]). We expect perturbative QCD calculations to 
yield predictions for the inclusive decay rates of B hadrons into charmonium states that are 
comparable in accuracy to the predictions for their semileptonic decay rates. 

Most previous calculations of the rate for charmonium production in B meson decays 
[§, ||, H II have been carried out under the assumption that the production mechanism 
is the decay at short distances of a b quark into a color-singlet cc pair plus other quarks 
and gluons, with the c and c having almost equal momenta and residing in the appropriate 
angular-momentum state. We will refer to this as the "color-singlet mechanism". It was 
assumed that the only nonperturbative input required in the calculation is the nonrelativistic 



cc wave function (for S states) or its derivative (for P states) at the origin. In this paper we 
point out that this assumption fails for P states, the failure being signaled by the presence 
of infrared divergences in the QCD radiative corrections. The divergences appear because 
there is a second production mechanism for P states that also contributes at leading order 
in a s , namely the decay at short distances of the b quark into a color-octet cc pair in an S 
state, plus other partons. We will call this the "color-octet mechanism". In addition to the 
derivative of the nonrelativistic wavefunction at the origin, calculations of P-wave production 
require a second nonperturbative input parameter, the probability for the color-octet cc pair 
to radiate a soft gluon and form a color-singlet P-wave bound state. The two production 
mechanisms are summarized below by factorization formulas which are accurate to leading 
order in v 2 , the square of the relative velocity of the charmed quark and antiquark, and to 
all orders in the strong coupling constant a s (M&). The hard subprocess rates appearing in 
these formulas are calculated to leading order in the QCD coupling constant. Estimates of 
the nonperturbative parameter associated with the color-octet mechanism are obtained from 
experimental data on the production rate of Xc states in B meson decays. We find that the 
color-octet mechanism may account for a significant fraction of the observed decay rate into 
Xd and that B mesons should decay into \c2 at a similar rate. 

The color-singlet and color-octet mechanisms for the production of P-wave quarko- 
nium states have analogs in the decays of these states. For S-wave resonances, the electro- 
magnetic and light-hadronic decays proceed through the annihilation at short distances of 
the heavy quark and antiquark in a color-singlet S-wave state. The decay rates are then 
proportional to the square of Rs(0), the nonrelativistic wavefunction at the origin, with co- 
efficients that can be calculated perturbatively as a series in a s (Mg), where Mq is the heavy 
quark mass. In the case of P-waves, one might expect that the decay should also proceed 
through the annihilation at short distances of the QQ pair in a color-singlet P-wave state. If 
this were the case, the decay rates would be calculable in terms of a single nonperturbative 
factor R' P (0), the derivative of the nonrelativistic wavefunction at the origin. While this 
expectation holds true for decay rates into two photons, explicit calculations of the light 
hadronic decay rates reveal infrared divergences ||. It has recently been shown that these 
infrared divergences can be systematically factored into a second nonperturbative parameter, 
which is proportional to the probability for the quark and antiquark to be in a color-octet 



S-wave state at the origin |7|]. These results can be summarized by rigorous factorization 
formulas for the decay rates of quarkonium states into light hadronic or electromagnetic final 
states ]7|, |[. The factorization formulas are valid to leading order in v 2 , where v is the typical 
velocity of the heavy quark, and to all orders in the QCD coupling constant a s (Mg). 

The inclusive production rates for charmonium states also obey factorization formulas. 
To leading order in v 2 , where v is the typical velocity of the charmed quark relative to the 
charmed antiquark, and to all orders in a( a (M&), the inclusive decay rate into ip satisfies the 
simple factorization formula that was assumed in previous work: |2], || |], |5|] 

T(b^ij + X) = Gi r\ (b -> ccfSO + X) , (1) 

where T\ is the hard subprocess rate for producing a color-singlet cc pair with equal momenta 
and in the appropriate angular- momentum state 23+1 Lj = 3 Si. It can be calculated per- 
turbatively as a series in the QCD running coupling constant ct s {M^). The nonperturbative 
parameter G\ in (|l|) is proportional to the probability for the cc pair to form a bound state 
and is related to the nonrelativistic wavefunction at the origin Rs(0): 

_ 3 \R S (0)\ 2 

A phenomenological determination that makes use of the electronic decay rate of the i/j gives 
G\ ~ 108 MeV. It can also be given a rigorous nonperturbative definition, so that it can be 
measured in lattice simulations of QCD ||. The production rate of the radial excitation if}' 
obeys the same factorization formula (0), but with a parameter G\ that is smaller by the 
ratio of the electronic decay rates of ip and ip''- G\ ~ 43 MeV. 

In previous work ||, the decays of b quarks into P-wave charmonium states were 
assumed to satisfy simple one-term factorization formulas like ([!]). Because of the color- 
octet production mechanism, such formulas are incomplete, and their breakdown is signaled 
by the appearance of infrared divergences in order a s . The correct factorization formulas for 
P-wave production rates have two terms: 

r (b — h c + X) = Ht f i (b -> ccCPi) +X,fi) +3 H'M f 8 (b -> ccCSq) + X) , (3) 



T(b^ X cj + X) = HiTi(b^cc( 3 Pj)+X,fx) + (2J + 1) H' 8 (fi) T 8 (b ^ cc( 3 S 1 ) + X) . 

(4) 

The factors I\ and r 8 are hard subprocess rates for the decay of the b quark into a cc pair in 
the appropriate angular-momentum state with vanishing relative momentum, the cc being 

in a color-singlet state for Ti and in a color-octet state for r 8 . They can be calculated 
perturbatively as series in at a (M&), with coefficients that are free of infrared divergences. 
Note that in (|4]) the only dependence on the total angular momentum quantum number 
J in the second term on the right side lies in the coefficient 2 J + 1. The nonperturbative 
parameters Hi and H 8 are proportional to the probabilities for a cc pair in a color-singlet 
P-wave and a color-octet S-wave state, respectively, to fragment into a color-singlet P-wave 
bound state. They can be given rigorous nonperturbative definitions || in terms of matrix 
elements in nonrelativistic QCD. The parameter Hi is directly related to the nonrelativistic 
wavefunction for the heavy quark and antiquark: 

^ 9 |^(0)! 2 (5) 

Its value has been determined phenomenologically from the light hadronic decay rates of 
the Xci and Xc2 to be Hi ~ 15 MeV 0. There is no rigorous perturbative expression for 
H' 8 in terms of the wavefunction Rp(r), since a cc pair in a color-octet S-wave state can 
make a transition to a color-singlet P-wave state through the radiation of a soft gluon, which 
is a nonperturbative process. The parameter H' 8 is not related in any simple way to the 
analogous parameter H 8 that appears in P-wave decays j7|, so it cannot be determined from 
data on the decays of P states. A phenomenological determination of H 8 can come only from 
a production process. Below we will use experimental data on \ c production in B meson 
decays to obtain a rough determination of H 8 . 

In the factorization formulas (^) and (|]), H' 8 and Ti depend on an arbitrary factor- 
ization scale fx in such a way that the complete decay rate is independent of fi. In order 
to avoid large logarithms of /x/Mj, in Ti, the factorization scale fx should be chosen to be 
on the order of M&. The scale dependence of H' 8 (fx) is governed by a renormalization group 



equation, which to leading order in a s (fi) is [|J 



The matrix element Hi is independent of the factorization scale \x. The solution to the 
renormalization group equation is therefore elementary. For example, for \i < M c , the 
solution is 



H' 8 (M b ) = H'M + 



16 h / a.(ji) \ + 16 ln fa s (M c 



27 fa \a s (M c ) 27 (3 4 \a s (M b 



Hi , (7) 



where (3 n = (33 — 2n)/6 is the first coefficient of the QCD beta function for n flavors of 
massless quarks. In the limit in which M b is very large, the contribution to H 8 from the 
perturbative evolution dominates, and one can estimate H' 8 (M b ) by setting a s (fi) ~ 1 and 
neglecting the constant H' 8 (n) in (0). Taking a s (Mb) w 0.20 and a s (M c ) w 0.31, we obtain 
H' 8 (M b ) ps 3 MeV. This should be regarded as only a rough estimate, since the physical value 
of M b is probably not large enough for the constant H' 8 (fi) to be negligible. 

The subprocesses in the factorization formulas (|3]) and (^) are decays of the b quark 
into cc plus other quarks and gluons. The dominant contributions involve the coupling 
of the b quark to ccs via an effective 4-quark weak interaction ||, which can, by Fierz 
rearrangement, be put into the form 

q ( 2C C 

C-weak = ~^=V cb V* s ( — ^ C7 M (1 - 7 5 )c ^(1 - 75)6 

+ (C + + C-) c 7m (1 - l,)T a c SY(l - l,)T a b^ , (8) 

where Gp is the Fermi constant and the V^'s are elements of the Kobayashi-Maskawa mixing 
matrix. The weak interaction that gives the Cabibbo-suppressed transition b — > ccd is 
obtained by replacing s by d and V* s by V* d in (]|). The coefficients C + and C_ in @ are 
Wilson coefficients that arise from evolving the effective 4-quark interaction mediated by the 
W boson from the scale M w down to the scale M b . To leading order in a s (M b ) and to all 



orders in ot s (M b ) ln(M^/Mb), they are 



C . W) M ) 12/23 . (10) 

Taking a s (M w ) = 0.116 and a s {M b ) = 0.20, we find that C + {M b ) « 0.87 and C_(M 6 ) w 
1.34. When a 6 quark decays through the interaction the first term produces a cc pair 
in a color-singlet state, while the second produces a cc pair in a color-octet state. The color- 
singlet coefficient 2C + — C_ is decreased dramatically by renormalization-group evolution, 
from 1 at the scale Mw, to 2C + — C_ ~ 0.40 at the scale M b . The color-octet coefficient 
(C+ + C_)/2 is increased slightly, from 1 at the scale Mw, to (C + + CJ)/2 »s 1.10 at the 
scale Mft. Since the dramatic suppression of 2C + — C_ is due to a cancellation between 2C+ 
and C_, it is sensitive to both the choice M& for the scale and to higher-order perturbative 
corrections to the Wilson coefficients. This sensitivity can be removed only by calculations 
beyond leading order. The Wilson coefficients C + and C_ have been calculated at next-to- 



leading order [ II| , but the calculations are meaningful only when combined with decay rates 



that are also calculated beyond leading order. 

It is convenient to express all the subprocess rates appearing in the factorization 
formulas ([l]), (§), and in terms of 



Gl /„ .Ml 



\V cb \ 2 —^MlM c I 1 - 4^ I . (11) 
1 cb| 144vr 6 c I M 2 ' 



To leading order in a s (M&), the color-singlet subprocess rates I\ can be extracted from 
previous calculations || H, !]• The sum of the two subprocess rates that contribute to the 
decay into the ip at leading order is 

f l (b ^ C c(3 Sl ) +s,d) = (2C + - C_) 2 (l + 8^§) f o , (12) 
where we have used |V^ S | 2 + \ V c d\ 2 ~ 1. The subprocess rate that contributes to the decay 



into Xci is 



fx (b - cc( i P 1 ) + s, d) = 2 (2C+ - C_) 2 (\ + 8^f) f . (13) 

The color-singlet subprocess rates that contribute to the production of h c , Xco, and Xc2 (the 
1 Pi, 3 Po, and 3 P2 states, respectively) vanish to leading order in a s because of the J PC 
quantum numbers of these charmonium states. 

The rate for direct production of ip in B meson decay has been calculated to next- 
to- leading order in a s (Ref. ||). Unfortunately, renormalization group effects were treated 
incorrectly in that calculation. The treatment of Ref. || was equivalent to using 2C + — 
C_ = (a s (M b )/a s (M w ))~ 2A / 23 and (C++C_)/2 = (a s (M b ) / a s {M w )f/ 23 for the color-singlet 
and color-octet coefficients at the scale M b . This treatment correctly reproduces the term 
proportional to a s ln(Mjy / 'M b ) in the order a s correction, but it fails to reproduce the leading 
logarithms at order a 2 s and higher. The results of Ref. || were presented only in graphical 
form, which prevents us from extracting the correct order-a s contribution to the subprocess 
rate for ip production. In using the leading order result (0), one should keep in mind that 
the next-to- leading correction proportional to a s {C + + CJ) 2 may be as important numerically 
as the leading term, which is proportional to (2C + — CJ) 2 . 

The color-octet subprocess rates T 8 appearing in (H) and (|) require a new calculation. 
The most straightforward method (although not the simplest) is to calculate the infrared 
divergent part of the rate for the decay b — > ccsg, with the cc in a color-singlet P-wave 
state after having emitted the soft gluon. This rate can be calculated in terms of Pp(0) by 



making use of a covariant formalism fl2|| . The infrared divergences come from the region 
of phase space in which the momentum of the radiated gluon is soft. In this region, the 
matrix element can be factored into the amplitude for b — > ccs, with the cc projected onto 
an S state, and a term that depends on the gluon momentum. The divergence comes from 
integrating over the phase space of the gluon. Imposing an infrared cutoff \x on the energy 
of the soft gluon, one finds that the divergence is proportional to ln(M&/ \i)H\. The final 



results for the infrared divergent terms in the decay rates are 



T(b^h c + sg) ~ - 2 \V cs \\C + + C^a s \n^ l -^lf , (14) 

T(b^ XcJ + sg) ~ (2J + l)A| Ks P(C + + C_) 2 a s ln^^|^ (l + 8^) f . 

(15) 

One can identify the corresponding infrared divergence in the perturbative expression for H' 8 
by neglecting the running of the coupling constant in ([]). The resulting expression for is 

H' 8 (M b ) ~ ^-a.ln^ F x . (16) 

Using this identification together with the expression (|5|) for Hi, we find that the color-octet 
subprocess rates defined in (^) and (|) are 

f 8 (b - cc( x 5o) + s, d) = I (C+ + C_) 2 f o , (17) 



f 8 (b - cc( 3 Si) + s, d) = I(C + + C_) 2 ^l + 8^ f 



(1* 



We now turn to the phenomenological applications of our results. Among the correc- 
tions to the factorization formulas for B hadron decays are the effects of spectator quarks 
and antiquarks, which are suppressed by powers of A /M& (Ref. p]). It is clear from decays 
of D mesons that spectator effects are much more important for total decay rates than for 
semileptonic decays. Our predictions for the inclusive decay rates into charmonium states 
should therefore be more reliable if they are normalized to the semileptonic decay rate, in- 
stead of the total decay rate. These branching ratios should be identical for B~ , B°, B s , 
and Aft, even if their lifetimes differ substantially. To leading order in a s , the semileptonic 



decay rate is |TT 



r(6 -> e~U e + X) = \V cb \ 2 ^-Ml F(M C /M b ) , (19) 



where F(x) = 1 — 8x 2 — 24x 4 lnx + 8x 6 — x 8 . Approximating the ratio M c /Mb of the heavy 
quark masses by the ratio Mr, /Mb ~ 0.35 of the corresponding meson masses, we find that 
the phase space factor is F(M C /M b ) rs 0.41 and that the semileptonic decay rate ( |H)D is 
|Kfe| 2 (M fe /5.3GeV) 5 (3.9 • 10~ n GeV). The Kobayashi-Maskawa factor \V cb \ 2 , as well as the 
extreme sensitivity to the value of the bottom quark mass M&, cancels in the ratio between 
the charmonium and semileptonic decay rates. 

The leading order QCD prediction for the inclusive decay rate into ip is obtained by 
inserting (|i"2"|) into (|1|). To take into account phase space restrictions as accurately as possible, 
we approximate the ratio M c /Mb of quark masses by M^ / (2Mb) ~ 0.29. Normalizing to the 
semileptonic decay rate, we obtain the ratio 

*W - m (t ~-- + i^ K M (2C+ - c - f W ■ < 20 > 

T(b — > e v e + A ) M 6 

Multiplying by the observed semileptonic branching fraction for B mesons of 10.7%, we find 
that the predicted inclusive branching fraction for ip is 0.23%. There are large theoretical 
uncertainties in this result because order-a s corrections to the color-singlet Wilson coefficient 
2C + — C- may be significant and because the perturbative corrections to the subprocess rate 
(|I2"D proportional to a s {C + + C_) 2 need not be small compared to the leading-order term, 
which is proportional to (2C + — CS) 2 . These uncertainties can be removed only by a complete 
calculation of the subprocess rate Ti(b — > ipX) to order a s . There are additional theoretical 
errors due to uncertainties in the quark masses M c and relativistic corrections of order v 2 , 
and corrections of order A 2 /M 2 . However these uncertainties are probably small compared 
to the perturbative errors. 

For the P states, we define ratios of the inclusive charmonium and semileptonic decay 
rates, analogous to the ratio R(ip) defined in (p0[). We take into account the phase space 
restrictions in the decay rate into charmonium as accurately as possible by using half the 
bound state mass for M c and the B meson mass for M&, so that M c /Mb varies from 0.32 to 
0.34 for the various P states. For the mass of the h c , we have taken the center of gravity of 



the XcJ states. The resulting leading order QCD predictions for the ratios are then 

R(h c ) * 14.7(C + + C_) 2 ^^, (21) 

R( Xc0 ) * 3.2 (C + + C-) 2 ^^ , (22) 

R( Xcl ) « 12.4 (2C+ - C_) 2 § + 9.3 (C + + C_) 2 ^^ , (23) 

ikZft Mb 

R( Xc2 ) » 15.3 (C + + C_) 2 ^^ . (24) 

These predictions are subject to the same uncertainties as the prediction for ip given in (p0|). 

The ratios predicted above are for the direct production of charmonium states in the 
decay of the B hadron. The branching fractions that are measured directly by experiment 
include indirect production from cascade decays of higher charmonium states. The branching 
fractions for direct production can only be obtained by making assumptions about the cas- 
cade decays. The inclusive branching fractions that have been measured are (1.12 ± 0.16)% 
for^ |TJ, (0.46 ±0.20)% for f |TJ, (0.54±0.21)% for Xcl [0, and an upper bound of 0.4% 
f° r Xc2 [0- The branching fractions for the cascade processes ip' — > tp + X, ip' — > Xc \ + 7, 
and Xci — ► V' + T are approximately 57%, 9%, and 27%, respectively JHJ]. If one assumes that 
there are no other cascade processes that give significant contributions, then the branching 
fractions for direct production by B meson decay are (0.71 ±0.20)% for ip, (0.46±0.20)% for 
ip f , and (0.50 ±0.21)% for Xc \. The ratio of the direct production rates of ip' and ip is consis- 



tent, within the large error bars, with the ratio of their electronic widths, which is 0.40 |13 



However the result for ip is several standard deviations larger than the prediction of (pop. 
The discrepancy could be due to order-a s corrections to the Wilson coefficient 2C + — C_ or 
to corrections to the subprocess rate (|12"D proportional to a s {C + + C ) 2 . Alternatively, it 
could be due to contributions from cascade decays into ip and ip' from higher charmonium 
states. Studies of the momentum distribution |1| and the polarization 0, [| of the ^'s could 
help determine whether cascades from states other than ip' and Xc \ are important. 



The measurement of the rate for B meson decay into Xc \ (Ref. [|T^]) makes it possible 
to extract a phenomenological value for the unknown matrix element H f 8 . Dividing the 



branching fraction for direct decay of the B meson into Xci by the semileptonic branching 
fraction (10.7 ± 0.5)% @, we find the branching ratio (J2|) to be 0.047 ± 0.020. With the 
choice Mi, = 5.3 GeV, the color-singlet term on the right side of (|23"D is 0.006. Attributing 
the remainder to the color-octet term, we find that ifg(Mb) « (4.8 ± 2.3) MeV. The quoted 
error is due to the experimental error in the branching ratio only, and does not include the 
theoretical error due to the potentially large perturbative corrections to the Wilson coefficient 

2C + — C- and to the color-singlet subprocess rate T\. One can also use ( p4|) together with 
the upper bound on the branching fraction into Xc2 to obtain the upper bound ifg(Mk) < 2.7 
MeV. This bound is consistent with the determination from Xci production, given the large 
error bar. The large theoretical uncertainty in our determination of H f 8 could be reduced 
dramatically by calculating the order-a s corrections to the P-wave subprocess rates. In the 
absence of these calculations, one should regard our determination of H' s as only a rough 
estimate. 

The color-octet mechanism has a dramatic effect on the pattern of the produc- 
tion rates for P-wave charmonium states. Taking the value if§(Mb) « 2.5 MeV, which 
is consistent with both the determination from decays into Xci and with the upper bound 
from decays into Xc2, we find that the relative production rates predicted by (|2l|)-(p4]) are 
h c '■ Xco '■ Xci '■ Xc2 = 1-3 : 0.3 : 1 : 1.3. A previous calculation of the decay rates of B mesons 
into charmonium states which did not take into account the color-octet mechanism, gave 
the relative rates h c : Xco '■ Xa '■ Xc2 = 0:0:1:0. The observation of Xc2 production 
at a rate comparable to that for Xci production would be a dramatic confirmation of the 
color-octet production mechanism and would provide a direct measurement of the matrix 
element H'%{M h ). 

P-wave production provides unique information on the production of heavy-quark 
pairs and on their binding into quarkonium. In this paper we have outlined a consistent 
factorization formalism for describing these processes in QCD. In addition to the familiar 
color-singlet production mechanism, the factorization formulas take into account the color- 
octet mechanism, in which a heavy quark- ant iquark pair is created at short distances in a 
color-octet S-wave state and subsequently fragments by a nonperturbative process into a 
color-singlet P-wave bound state. The inclusive production rates for all four P states can 



be calculated in terms of two nonperturbative inputs: a parameter H\ that is related to 
the derivative of the wavefunction at the origin and a second parameter H' 8 that gives the 
probability for the fragmentation from the color-octet S-wave state. We have applied this 
formalism to the production of P-wave charmonium in decays of B hadrons and have used 
experimental data to obtain a crude estimate of the parameter H' s . An accurate determina- 
tion of Hg requires perturbative calculations beyond leading order as well as more accurate 
experimental data on \c production in B meson decays. The accurate determination of 
this parameter would be very useful because the production rates for P-wave charmonium 
states in photoproduction, leptoproduction, hadron collisions, and other high energy physics 
processes satisfy factorization formulas involving the same two nonperturbative parameters 
Hi and H' 8 that appear in B meson decay. Knowledge of these two parameters would, in 
principle, allow one to compute all of the inclusive P-wave charmonium production rates. 
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